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ABSTRACT 

We have computed a series of realistic and self-consistent models of the emitted spec- 
tra of Hn galaxies. Our models combine different codes of chemical evolution, evolutionary 
population synthesis and photoionization. The emitted spectrum of Hn galaxies is reproduced 
by means of the photoionization code CLOUDY, using as ionizing spectrum the spectral en- 
ergy distribution of the modelled Hn galaxy, which in turn is calculated according to a Star 
Formation History (SFH) and a metallicity evolution given by a chemical evolution model 
that follows the abundances of 15 different elements. The contribution of emission lines to the 
broad-band colours is explicitly taken into account. 

The results of our code are compared with photometric and spectroscopic data of Hn 
galaxies. Our technique reproduces observed diagnostic diagrams, abundances, equivalent 
width-colour and equivalent width-metallicity relations for local Hn galaxies. 
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1 INTRODUCTION 

Hn galaxies, a subset of Blue Compact Galaxies (BCG), are gas 
rich dwarf galaxies whose optical spectra are dominated by strong 
and narrow emission lines. They are currently experiencing in- 
tense star formation in small volumes. The emission lines are pro- 
duced by gas ionized by the young massive stars. Observations 
indicate that Hn galaxies are in general metal poor systems. This 
fact in addition to their very young stellar populations has been 
known for some time and has led to the proposal that these sys- 
tems are very young, su ffering their first burst of star formation 
jSargent&Searielll97(3h . 

It seems, however, that these galaxies are not as young as it 
was originally thought, since recent photometric observations indi- 
cate that they host stellar populations which a re at least 10 7 - 10 s 
a. old, reaching in some cases a few Ga I Telles et"aH 1 997 ; 



Legrandl. [200f3;lTolstoyl. l2003l : ICairos et all 120031 : Ivan Zee et a! 



2004lThuan & IzotovmOOsTT There is now wide agreement about 



the existence of underlying populations with ages of several Gas 
in BCGs and/or Hu galaxies. Even in the most metal-poor galaxy 
known, Izwl8, the best candidate to primordial .g aJaxy, there is an 
under lying population of at least some 10 8 -10 9 a l llzotov & Thuari 
120041) . 
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However there is no consensus about their history of star for- 
mation. Two scenarios are postulated: 

a) The star formation takes place in short but intense episodes 
separated by long quie scent periods of null or low activity 
dBradamante et alill998h . 

b) The star formation is continuous and of lo w intensity alon g 
the galaxy life, with superimposed sporadic bursts (Legrand. 2000). 

In both scenarios, a stellar population of intermediate age 
contributes to the observed luminosity. Indeed, recent work takes 
into account the existence of an underlying stellar population 
and combine photo-ionization codes with evolutionary synthe- 
sis models ( hereinafter ESM). Ex a mples of this procedure ca n 
be found in iGarnett et al.l i ll 9991). iMas-Hesse & Kuntr] i 19991) . 
IStasinska & Izotovl J2003h and iMov et alj d200ll) . In general the 
ESM so far do not include the chemical evolution. Given the low 
metal content of Hu galaxies the metallicity effect can be of con- 
siderable relevance. 

iTerlevich et al. I J2004l) have shown that the equivalent width 
of the H0 line, EW(H/3), which can be taken as an age indicator 
for the ionizing population, and the (U-V) colour in Hu galaxies 
show an anti-correlation which cannot be interpreted in terms of 
a single stellar population (SSP). At a given EW(H6) the data are 
displaced to colours redder than those predicted by the models. The 
most likely explanation is that the ionizing population, that is the 
most recent burst of star formation, is overimposed on an older stel- 
lar population evolved enough as to produce a redder (U-V) colour. 
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There is also an inverse correlation between EW(H/?) and the gas 
oxygen abundance, 12 + log(0/H), which could be interpreted in 
terms of galactic chemical evolution: as the average age and the to- 
tal mass of the stellar population increases , the EW(H/3) decreases 
and the mean metallicity o f the interstellar medium increases. 

A number of works dShiet all 120061 ,and references therein) 
have computed purely chemical evolution models (hereinafter 
CEM) for the study of BCG and dwarf irregular galaxies. Most 
of them assume that the star formation occurs in bursts and in- 
clude the effects of galactic winds and/or gas infall. However most 
of them limit t he study to the evolution of nitrogen and oxy- 
gen abundances faenrv et alj|2000l: lLarsen et al.ll200 1|) and /or the 
luminosity-me t allicit y relation~ ]lVIouhcine & Contin 3 120021) . Only 
IVazquez et alj J2003h used the informat i on com ing from the chem- 
ical evolution models from lCarigi eTai] d 2002h to perform the next 
step and combine chemical and spectral evolution. Their models 
exclude however the early stages of evolution, i.e. during the nebu- 
lar stage when the most massive stars dominate the energy output. 

A code combining chemical evolution, evolutionary synthesis 
and photo-ionization models has not yet been developed and used 
for the spectral analysis of BCG and Hn galaxies. 

In this work, we include the chemical evolution model results 
in the computation of the spectral energy distributions which are 
used as ionizing source for a photo-ionization code. 

Our method allows the simultaneous use of the whole avail- 
able information for the galaxy sample concerning on the one hand 
the ionized gas - emission lines intensities and equivalent widths, 
elemental abundances, gas densities - which defines the present 
time state of the galaxy, and on the other hand spectro-photometric 
parameters - colours, absorption spectral indices — defined by the 
stellar populations and their evolution with time and thus related to 
the galaxy star formation history. This is done in a self-consistent 
way, that is using the same assumptions regarding stellar evolution, 
model stellar atmospheres and nucleosynthesis, and using a realis- 
tic age-metallicity relation. 

The basic codes employed are described in Section 2. Section 
3 shows the corresponding results of each kind of model, and the 
whole combined results about galaxy evolution. These results are 
discussed in Section 4. Finally, our conclusions are summarized in 
Section 5. 



2 MODEL DESCRIPTION 
2.1 Chemical evolution 

The code used for this work is a generalization o f the code used for 
the Solar Neighborho od bv|Ferrini et al.l 1 1992h and applied to the 
whole Galaxy disc in iFerrini et al. 1 1994 ). That model described 
the galaxy as a two-zone system (halo and disc) in which the disc 
is a secondary structure formed by the gravitational accumulation 
of gas from the halo. In the present work, however, all the gas is 
assumed to be within an only region from t = 0, that is, the collapse 
time scale has been eliminated as an input parameter of the code. 
The matter is assumed to be in different phases: 



- A stellar population, where we distinguish the stars able to 
create and eject heavy elements to the interstellar medium (M > 
4M Q ), and the stars which only eject hydrogen and He (M < 4M Q ). 

- Stellar remnants, as the endpoint of star evolution, that act as 
a matter sink, removing mass from chemical evolution. 

- Interstellar diffuse gas out of which stars are forming follow- 



ing a simple Schmidt law. No molecular component is explicitly 
considered in the present work. 

A bursting star formation has been assumed taking place in 
successive bursts along the time evolution. At any time, the avail- 
able gas for each burst is the sum of the gas left after the previous 
burst of star formation and the gas ejected by massive stars. 

According to this description, the time dependence of the total 
mass fraction in each phase and the chemical abundances in the 
ISM are determined by the interactions between these phases and, 
as a consequence of the evolution, we get the star formation rate. 
The basic equations to study the behavior of the mass of stars and 
gas are, therefore: 



dM 
~~dt 

dM s 

~ 

dM g 

~dT 
M 





-*P + E 

M s + M e 



(1) 

(2) 

(3) 
(4) 



where M s (t) is the mass in stars, M g (t) is the mass of gas, M is the 
total mass of the system, v P(t) the star formation rate and E(t) the 
ejection rate of mass from the stars to the ISM. A part of the resti- 
tuted gas consists of enriched material, E z (t), which is created in 
the interior of the stars and is ejected when they die. To compute 
both quantities, ejected gas and ele ment production, we us e the 
Qij matrix formalism introduced by iTalbot & ArnetJ d 19731) and 
widely used in classical chemical evolution models. It is very use- 
ful to treat elemental abundances which increase at different rate 
in the interstellar med ium, that is, the non-solar elemental ratios 
jPortinari et alll 19981). In the presen t work we adopt the same ma- 
trix prescription as in IGavilan et alj l l2006l) . whic h is an updated 
versio n from the one given in lFerrini et al .1 i ll 9921) and iGalli et al 



19951) calculated following the method described bv lPortinari et al 



In order to calculate these matrices it is necessary to consider 
the nucleosynthesis of stars, their initial mass function (IMF) and 
their end either as type I or II supernovae, or through quiet evo- 
lution. The mean stellar lifetimes, final states, and nucleosynthesis 
products, depend on the stellar mass which makes the IMF to play 
an important role in the chemical evolution model. The a dopted 
IMF is assumed to be universal in spac e and constant i n time 1 Wvsel 
ll997l:IScalolll998l : lMever et alll2000h and taken fromlFerrini et al l 
Jl990h . This IMF is very similar to Scalo's law dScalolll986l) and 
in g ood agreement w ith the recent expressions from lKroupal ( l200lh 
and lChabrieJ ( l2003T). Nucleosynthes i s yield s for massive stars have 
been taken from lWooslev & Weaver! d 19951) . For low and intermedi- 
ate mass stars, we have used the yields from IGavilan et all d2005l) . 
The combination of these sets of stellar yields, with our assumed 
IMF producing the required yields per stellar genera tion, has been 
compl e tely revised and calib rated with the MWG in IGavilan et alj 
d2005l) : IGavilan et alj J200r3) . where other stellar yields sets have 
been used and analyzed t oo. It also reproduces o bservations of spi- 
ral and irregular galaxies dMolla & Diazll2005l) , obtaining reason- 
able results even fo r nitrogen abundances in low-metallicity objects 
dMolla et al .. 2006). For type I supernova ex plosion releases , mode l 
W7 from lNomoto et al]dl984l) , as revised bv liwamoto et alld 19991) . 
has been taken, while the rate of this type of SNIa explosions are 
included through a table give n by Ruiz-Lapuente (private commu- 
nication), following works o f'R uiz-Lapuente et a il booot) . 



Modelling the evolution ofHu galaxies 3 



We have assumed that the gas is used to form stars with a given 
efficiency called H: *F = H ■ M g . If there would not be gas ejection 
from massive stars, as it would be the case for instance in the firsts 
moments of the evolution, (E = 0), then: 



dt 
dt 



HM g 
-HM„ 



(5) 
(6) 



That is, the consumed gas rate would be a function that de- 
pends on the available gas at a given time, Q being a decreasing 
function of time through the parameter H, which defines the star 
formation efficiency. Then: 



dM. 

U7 = - Hdt ^ M ^ 



M g0 e 



(7) 



A change in the star formation rate implies a change in the 
value of the parameter H or efficiency: 



In 



H 



At 



(8) 



We have run models considering the star formation as a set 
of successive bursts in a region with a total mass of gas of 100-10 6 
M and a radius of 500 pc (see section 3.2 for details). Every galaxy 
experiences 1 1 star formation bursts along its evolution of 13.2 Ga, 
that is one every 1.3 Ga. In each burst a certain amount of gas is 
consumed to form stars. We have computed two types of models. 
In the first one, which we call burst models (models B), the same 
efficiency H is assumed for every burst while in the second one 
(models A) we assume an attenuated bursting star formation mode , 
i. e. a variable H is assumed along the time evolution. For each type 
we have computed two models with different initial star formation 
efficiency: 

- model 1 , that uses ~ 64 per cent of the gas for star formation 
in each time step. Thus, following Eq.(8), 1/H = 0.5 Ma 

- model 2, that uses ~ 33 per cent of the gas to form stars in each 
time step, which implies 1/H = 1.25 Ma 

These percentages are maintained in models B, which means that 
every burst will transform the corresponding percentage of the 
available gas (taking into account the gas ejected by massive stars) 
in stars. In Models A, these values correspond to the initial efficien- 
cies and are reduced in the successive bursts by a factor n, where n 
is the number of the burst. 

The code solves the chemical evolution equations to obtain, in 
each time step, the abundances of 15 elements: H, D, 3 He, 4 He, C, 
13 C, O, N, Ne, Mg, Si, S, Ca, Fe, and nr where nr are the isotopes 
of the neutron rich elements, synthesized from 12 C, 13 C, 14 N and 
16 inside the CO core. We have taken a time step, At= 0.5 Ma0, 
from the initial time, r = 0, up to the final one, t = 13.2 Ga. Also, at 
each time step, the star formation rate and the mass in each phase 
- low mass, massive stars and remnants, total mass in stars created, 
and mass of gas- are computed. This amounts to a total of 26348 
stellar generations or SSPs, each one with its corresponding age 
and metallicity. 



2.2 Evolutionary synthesis 

Once the chemical evolution models have been calculated, the spec- 
trophotometric properties can also be obtained. In th is work we 
have used the code from IGarcfa- Vargas et al.l dl995dlbl [ 1998) as 
updated in lMolla & Garcia- Vargas! d2000f) . The required inputs are 
the stellar isochrones and the stellar atmosphere models for the 
individual stars. In that work we used the s et of 50 isochrones 
from the Padova group from iBertelli et al. I dl994l) for ages be- 
tween 4 Ma and 20 Ga, for 5 metallicities between 1/50 and 
2.5 Z (Z=0.0004, 0.004, 0.008, 0.02 and 0.05). To this set of 
isochr ones we have added othe r 18 isochrones from lBressan et all 
dl993h andlFagotto et al.Ul994h with ages from 0.5 to 4 Ma used in 
IGarcfa- Vargas et alT i 1995al) and lGarcia- Vargas et al. I dl998h for the 
youngest stellar populations of 4 metallicities Z=0.004, 0.008, 0.02 
and 0.05. That means that for the lowest metallicity, Z = 0.0004, 
only isochrones older than 4 Ma are available. 

In both cases, isochrones give the number of stars in each 
phase assuming that stellar clusters were formed with a standard 



Salpeter IMF with nti ow = 0.6 M Q and m„ 



100 M . This IMF, 



although using the same lower and upper limits, is not exactly the 
one used in the multiphase chemical evolution code, what implies 
some differences in the resulting spectra, mainly for the oldest and 
the youngest stellar populations. They are, however, smaller than 
10 per cent (Garcia- Vargas et al. 2008, in preparation), and there- 
fore we assume that this inconsistency will not produce very im- 
portant differences in the final results. 

The emergent spectral energy distribution (SED) was syn- 
thesized by calculating the number of stars in each point of the 
H-R diagram and assigning to it the most adequate s t ellar a tmo- 
sphere model . The models of IClegg & Middlemassl Jl987l) and 
iLeieune et al.l i 19971) have been used for stars with T e ff ^ 50000 
K (last evolutionary stages of massive stars) and 2500 ^ T e ff < 
50000 K, respectively. 

To each SED computed for a stellar cluster, it was added 
a continuum nebu l ar em ission contribution as explained in 
IGarcfa- Vargas et alj Jl998l) . The gas is assumed to have an elec- 
tron temperature, T e which depends on Z. The assumed values are 
T f =18000, 11000, 9000, 6500 and 4000 K for Z=0.0004, 0.004, 
0.008, 0.02 and 0.05, respectively, selected accord ing to the aver- 
age value obtained by IGarcfa- Vargas et al .1 l ll995bl) . The free-free, 
free-bound emission by hydrogen and neutral helium, and the two 
photon hydrogen -continuum h ave been calculated by means of the 
atomic data from lAlleJ J 19841) and lFerland Jl98Cl) . 

The SEDs for the complete sets of 68 (50 for Z=0.0004) 
isochrones were calc u lated with the same code as in 
iMolla & Garcia- Vargasl J2000h for the cited metallicities. 
Therefore we finally have a set of 68 SEDs for SSP of different 
ages for 4 metallicities each, and 50 for ages older than 4 Ma for 
Z=0.0004, which we use as a spectral library. 

Once the SEDs for SSPs are computed, they must be con- 
volved with the star formation history (SFH) in order to calculate 
the corresponding SED for a region where more than a burst take 
place: 



L A (t)= f S A (r,Z(t'W(t')dt' 
Jo 



(9) 



where r = / - /' is the age of the stellar population created in a 



When E + this expression is modified, as we will show in Section 3. 
2 This small time step allows to take into account the fast evolutionary 
phases of the most massive stars considered. 



3 Although in that work, applied to old elliptical galaxies, the SEDs cor- 
responding to ages younger than 4 Ma and those for Z=0.0004, were not 
used 
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time ?' and 5,! being the SED for each SSP of age r and metallicity 
Z reached in that time f . A SED from the SSP library, S A , must 
assigned to each time step according to its corresponding age and 
metallicity taking into account the SFH, ¥(0, and the AMR, Z(f), 
obtained from the chemical evolution model, to finally calculate 
L,i(0 by the above integration. However the metallicity changes 
continuously while the available SEDs have only 4 or 5 possible 
values. We have, therefore, interpolated logarithmically between 
the two SSP of the same age r closest in metallicities to Z(t') to 
obtain the corresponding S a (t, Z(t')). The final result is the total 
luminosity at each wavelength A. 

When this process is applied to the stellar population created 
by the first burst two problem arise. On the one hand the initial 
metallicity is Z= and, although during the first Ma after the cre- 
ation of the first stars the metallicity increases, it does not reach the 
minimum Z=0.0004, and therefore we must extrapolate with the 
available SEDs of the youngest SSPs. Due to the uncertainties in 
the evolution of the quasi-zero-metallicity stars, we have preferred 
to assign to this first burst metallicity a minimum value of 0.0028, 
smaller than the one reached at the end of the first burst but still 
valid for extrapolating without creating instability numerical prob- 
lems. On the other hand we must to use two different pairs of metal- 
licities to extrapolate with the SEDs of SSPs: Z=0.004 and 0.008 
before 4 Ma and Z=0.0004 and 0.004 after 4 Ma, since SEDs for 
the most-metal-poor and youngest SSPs are not available. We have 
checked that this method does not produce discontinuities and that 
the resulting spectra behave smoothly. In this way we obtain the 
SED corresponding to the whole stellar population, including the 
ionizing continuum proceeding from the last formed stellar popu- 
lation. 



2.3 Photo-ionization 

The photo-ionizatio n code CLOUDY, in its 96.0 version 
dFerland et all 1 19981) , has been used in order to obtain the emission 
lines. The gas is assumed to be ionized by the massive stars be- 
longing to the current burst of star formation whose SED has been 
previously calculated by the combination of the chemical and evo- 
lutionary synthesis models. The shape of the ionizing continuum is 
defined by the pair of values (y(Ryd), logvL v ) (Eq.[9]l. 

We have assumed the emitting gas to be located at a distance 
R = 500 pc at the begi nning of the evoluti on. This size is charac- 
teristic of Hn galaxies dTelles et"all I 199% . This radius is subse- 
quently adjusted as to keep a constant stellar density configuration. 
This results in all cases in a plane-parallel geometry. 

A given model is characterized by its ionization parameter U 
defined as: 

V = Q(H)l4ncn H R 2 (10) 

where n e is the electron density, assumed constant for simplicity 
and equal to 100 cm 4 , and Q{H) is the number of Lyman ionizing 
photons, calculated from the cluster SED, which reach the gas at a 
velocity c. Finally, the gas elemental abundances are those reached 
at the end of the starburst previous to the current one, as calculated 
from the chemical evolution model, except for Na, Ar and Ni which 
are not computed in the model and are scaled to the solar ratio 
dAsplund et alll2005h . 



3 RESULTS 

3.1 Chemical evolution 

3.1.1 Star formation rates 

The star formation rate is one of the main results of the different 
models since it drives the behaviour of all the other quantities. In 
Fig.[T]the star formation history for each of the 4 computed models 
is shown. In all of them the first burst is strong, while the succes- 
sive bursts are less intense since the amount of gas available has 
decreased in spite of the ejected gas produced by the massive stars. 
In models A (right panels), the first burst has the same star forma- 
tion rate as in Models B (left panels), since the initial efficiency 
is the same for both. The successive star formation episodes have a 
smaller star formation rate in models A than the corresponding ones 
in m odels B due to the a ttenuation factor included in the inputs. 

iHovos et alj j2004h have estimated the star formation rates for 
the sample of 39 lo cal Hn galaxies of different morphologies from 
iTelles et all d 1997b. which were selected from the most luminous 
of the lTerlevich et alj dl99lh catalog and with H/3 equivalent widths 
between 30 and 280 A . These star formation rates were obtained 
from the Ha luminosity, by assuming T e = 10 4 K and case B re- 
combination. They are in the range from 0.5 to 7 MgCT 1 , marked 
by dotted lines in Fig. Q] indicating that these Hn galaxies have 
high star f ormation rates. The se findings are in agreement with 
those from iTavlor et al 1 d 19941) who found star formation rates in 
the range of 10~ 4 to 10 MqcT 1 . Our models produce SFR between 
the observed values after the sixth burst, which corresponds to 6.5 
Ga after the beginning of star formation, except model B2, where 
every burst before the ninth-tenth show higher SFR than observed. 

Since the star formation consumes gas, the averaged gas mass 
decreases with time following Eq.(6), as can be seen in Fig. [2] The 
chemical evolution code, however, takes into account the mass of 
gas ejected from massive stars during their evolution and, therefore, 
the gas mass increases during a given burst. On the other hand, the 
mass of stars which increases in each burst and globally due to 
the star formation, decreases during each burst as the most massive 
stars end their lives. It is interesting to note that models with smaller 
efficiencies, H, in the lower panels of Figs. Q] and [2] consume less 
gas in each burst and, as a consequence, the star formation rate re- 
mains higher, and declines more slowly, than in models with higher 
values of H shown in the upper panels. 

3.1.2 Elemental abundances 

The element content of the gas is usually measured through the 
oxygen abundance, given as 12 + log(0/H). Oxygen contributes 
more than 40 per cent to the mass in metals and is mainly created 
within massive stars. Since these stars are very short lived, their 
products are ejected very rapidly to the interstellar medium (ISM) 
and hence the value of its oxygen abundance increases with time, as 
shown in Fig. [3] In that figure we see how O/H increases abruptly 
when a burst takes places remaining more or less constant between 

each two of them. 

Observati o nal d ata from iTerlevich et al. dl99lh and 
IHovos & Diazl d2006l) show that most of the Hn galaxies 
have metallicity distributions between 7.5 < 12 + log(0/H) < 8.5. 
These limits are shown as dashed lines in Fig. [3] The first limit is 
reached with the first burst in all cases. It can be seen that models 
of type 1 (dotted lines) show oxygen abundances higher than the 
observed upper limit. The efficiency of the first burst is so high that 
a large number of massive stars are formed, and hence the oxygen 
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Figure 1. Star formation history for: a) burst model 1 (Bl); b) attenuated burst model 1 (Al); c) bur st model 2 (B2); an d d) attenuated burst model 2 (A2). 
Dotted lines represent the upper and lower limits to the SFR estimated for BCG and/or Hii galaxies bv lHovos etall 12004 . 




Figure 2. Time evolution of the gas and star mass for models: a) Bl; b) Al; c) B2; and d) A2. 



abundance reaches rather high values in a very short time. Also 
the oxygen abundance of model B2 reaches values higher than 
shown by data from the third burst onwards (12 + log(0/H) > 8.5). 
Only model A2 shows oxygen abundances within the range of Hn 
galaxies during the whole evolution since the attenuation of the 
bursts keeps the star formation lower than in the other three cases. 
The star formation efficiency taken for model A2 results to be 
an upper limit for HII galaxies in oxygen abundances and it may 
reproduce the most metal rich systems. 

Abundance ratios also give information about the star forma- 
tion process, essentially the time scale or duration of the bursts, 
when the production of the elements involved takes place in stars 
in different mass ranges. This is the case for the N/O ratio. Oxygen 
is very quickly ejected by the most massive stars. Nitrogen, how- 



ever, may be created in stars of all masses. Moreover, the produc- 
tion of nitrogen may take place in two ways: 1) the main process 
requires another element like carbon or oxygen, to synthesize N 
through the CNO cycle, thus N shows, at least in part, a secondary 
behavior. In that case the N abundance has to be proportional to 
the initial abundance of heavy elements and hence the N/O ratio 
grows with metallicity. There is, however, a proportion of Nitro- 
gen w hich must have a primary origin, as suggested by the data in 
Fig.l4l dlzotov & Thuanl . [T99^ ; llzotov eTail l2005t [2000) shown by 
(green) dots which exhibit a constant N/O ra tio. Nitrogen may be 
produced by low and intermediate mass stars jGavilan et alll200r3) 
and therefore its contribution may appear in the ISM a certain time 
after the massive stars have died. In fact, N grows more slowly than 
the oxygen abundance in the range 7.6 < 12 + log(0/H) < 8.2, 
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Figure 3. Time evolution of the oxygen abundance, given as 12 + log(0/H), for the 4 kinds of computed models: a) models of type B burst modes) and b) 
models of type A attenuated mod els) with dotted lines c orre sponding to models of type 1 and solid lines corresponding to models of type 2. The dashed lines 
limit the range of data taken from lTerlevich et alj jl99ll) and lHovos & Di'azl j2006h . 



which may be explained by this primary behavior (see M oTla et all 
l2006t for details). In Fig.H we can see the oscillating behavior of 
the N/O ratio shown by our computed models due to the successive 
bursts of star formation followed by quiescent periods. At the be- 
ginning of a given burst, the N/O ratio decreases, as the O is created 
and ejected by the most massive stars in the burst; then it increases 
during the quiescent periods between bursts as the low and inter- 
mediate mass stars die and eject their synthesized nitrogen while 
the oxygen abundance remains constant. 

According to the results of this section, model A2 is the one 
better reproducing the observations, since the total mass, the gas 
fraction, the star formation rate and the oxygen abundances are 
within the range shown by data. Therefore, from now onwards, in 
what regards the ionizing continuum and emission line properties 
only the A2 model will be analyzed. 



3.2 Evolutionary stellar population synthesis 

The convolution of the results of the chemical evolutionary code 
with the single stellar populations of the synthesis evolutionary 
code gives us directly the SED, h A , in time steps of logt= 0.05. 

From this SED the number of hydrogen ionizing photons, 
coming from hot and young stars of the first spectral types (O-B), 
can be calculated as: 



Q(H) = 




where L v is the continuum luminosity at frequency v. 

Fig. [5] shows the number of hydrogen ionizing photons as a 
function of time for model A2. In the first burst, with a total mass 
of ~ 30 ■ 10 6 M G forming stars, the number of massive ioniz- 
ing stars is high and so it is the number of ionizing photons. The 
number of ionizing photons per unit stellar mass in our models is 



log(2(H)/M ) ~ 46.8 at the starting of the first burst, correspond- 
ing to a Zero Age Stellar Population (ZASP) of very low metallic- 
ity and then decreases as the cluster ages. This initial number of 
ionizing photons per unit stellar mass is almost the same for the 
successive bursts although, decreasing slightly as the abundance 
increases. 

Within each burst the number of ionizing photons decreases 
by almost two orders of magnitude from 0.5 to 10 Ma, as show in 
Fig. [6] and by almost 8 orders of magnitude by 100 Ma. Essentially 
no ionizing photons are available during the quiescent inter-burst 
periods. At the beginning of each successive burst, the number of 
ionizing photons raises abruptly but reaching a value progressively 
lower as less gas is available to form stars. Also, as the mean gas 
metallicity increases with ti me, the number of ionizing ph otons per 
unit stellar mass decreases dGarcia- Vargas et ailll995aB bl). This ef- 
fect increases with the age of the ionizing cluster reaching a factor 
of about 5 at 10 Ma. 

The average number of io nizing photons observ ed for local Hn 
galaxies is log (£>(#)) = 53.3 dHovos & Drazll2006l) . This number 
of ionizing photons, marked as a dotted line in Fig. [5] is reached by 
model A2 at around the fifth or sixth burst of star formation, that is 
for the burst occurring 5.2-6.5 Ga after the beginning of the forma- 
tion of the galaxy, when the predicted star formation rate also fits 
the data (see Fig. [T). On the other hand, the average derived ion- 
ization parameter for the local Hn galaxy sample is logU = -2.5 
which in our models corresponds to an age of the ionizing popula- 
tion of about 6.5 Ma. It should be noted that during the first 7 Ma of 
a burst the number of ionizing photons remain essentially constant 
(see Fig. [6J. 

Therefore, if the current burst is identified with redshift z = 0, 
the formation of the galaxy would have taken place at a r edshift 
z ~ 0.7 ( using H Q = 74 km.s'.Mpc -1 , iMacri et al.Ll2006l) . Also 
the observed number of ionizing photons for LCBGs at a redshift 
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Figure 4. The logarithmic N/O ratio as a function of the oxygen abundance, given as 12 + log(0/H), for our 4 kinds of computed models: a) type B models; 
b) type A models. The (black) dashed lines limit th e range of data shown by extragalacti c Hn regions . The (green) dots and ( magenta) squares c o rrespond to 
the de rived values for low metallicity galaxies from llzotov & Thuanl <1999l) ; llzotov et al.l {2005, 2006); IHovos & Di'aJ fe006t) ; lLiang et alj l2O06l) ; lNava et alj 
120061) . while the large triangle is the value for the low est metallicity gala xy in the sample SBS 0335-052W. The long dashed (cyan) lines represents the 
averaged trend shown by chemical evolution models from Molla et al. (2006) for spiral and irregular galaxies with continuous star formation. The arrows mark 
the solar values. 
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Figure 5 . The time evolution o f the hydrogen ionizing photons, Q(H), for model A2. The dotted line marks the mode of the distribution for local Hii galaxies 
found bv lHovos & DiaJ <2006l) . Each coloured line indicates the evolution within a burst. Different colour lines are used for different bursts. 



0.46 < z < 0.7, estimated as log(Q(H)) ~ 54.6 l lHammer etall 
1200 lb . would correspond to the value of the first burst of our mod- 
els. According to these results, our model seems consistent with ob- 
servations at these intermediate redshifts reproducing at the same 
time those corresponding to the more luminous local Hn galaxies. 
The question of if these galaxies are young galaxies, experi- 



encing their first burst of star formation or if, on the contrary, there 
exists an underlying population formed some Gas ago, may be ad- 
dressed by multi-colour photometry. Although a direct detection of 
the underlying population in LCBG or Hn galaxies is complicated 
because the observed light is dominated by massive stars, the ef- 
fect of an underlying stellar population should be easily seen in 
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Figure 6. The time evolution of the number of ionizing photons for each of the successive bursts during the 1.3 Ga inter-burst period. The burst number 
increases downwards. Different colour lines are used for the different bursts. 



the observed colours of these galaxies, since any star forming burst 
previous to the currently observed one will contribute substantially 
to the total continuum luminosity at the different wavebands. 

In Fig. [7] we represent the time evolution of the continuum 
U-B, B-V, V-R and V-I colours for our model A2 along the first 
10 Ma after each burst. The first burst -solid black thick line- fol- 
lows the expected evolution of a single stellar population of very 
low metallicity, remaining bluer than during the subsequent bursts 
where the effect of the accumulated continua from the previous star 
formation episodes makes colours become redder than expected for 
a SSP. Even though they are also starbursts in themselves, their star 
formation efficiency is much lower than that of the first burst, and 
so is their contribution to the total continuum luminosity which re- 
sults in redder colours. Furthermore, the metallicity, almost zero for 
the first burst, increases up to Z = Z Q /5 already during the second 
burst which also has an reddening effect over the colours. 

3.3 Photo-ionization models 

The calculated SED in each time step can be used as input ion- 
izing sources in the photo-ionization code to predict the emission 
lines intensities which provide information about the youngest stel- 
lar population which dominates the final spectra. 

Fig.[8]shows the relation between the [OIII]AA 4959,5007/H/3 
and [OII]/i/i 3727,3729/HjS ratios, which constitutes an excitation 
diagnostics for ionized nebulae. In this diagram, the most excited 
objects lie up and to the left while the objects with the lowest exci- 
tation are at the bottom right. 

The data shown in the graph have been ex tracted mainly 
from t wo main sources. First, the compilation from lHovos & Dfaa 
d2006h which provides emission line measurements, corrected for 
extinction, published for local Hn galaxies. The sample comprises 
450 objects and constitutes a large sample of local Hn galaxies with 
good-quality spectroscopic data. The sample is rather inhomoge- 
neous in nature, since the data proceed from different instrumen- 
tal setups, observing conditions and reduction procedures, but have 
been analysed in a uniform way. Data for these sample objects in- 
clude the emission line intensities of: [OII]^ 3727,29 A , [OUTJ/U 
4959,5007 A , and [NII]/U 6548,84 A , all of them relative to H/3, 
and the equivalent widths of the [Oil] and [OIII] emission lines - 
EWQOII]) EWQOIII]) -, and the H/3 line, EW(H£). As a second 



source, we have used the metal poor galaxy da ta from the Data Re - 
lease 3 of Sloan Digital Sky Su rvey, taken from Izotov et alj J200^> . 
The Sloan Digital Sky Survey l lYork et all[200oF constitutes a large 
data base of galaxies with well defined selection criteria and ob- 
served in a homogeneous way. The SDSS DR3 lAbazaiian et al .1 
d2005T) provides spectra in the wavelength r ange from 3800 to 9 300 
A for ~ 530000 galaxies, quasars and stars. Ilzotov et al .1 fcOOtj) ex- 
tracted ~ 2700 spectra of non-active galaxies with the [OIII]/i4363 
A emission detected above lcr level. This initial sample was fur- 
ther restricted to the objects with an observed flux in the H/3 emis- 
sion line larger than 10~ 14 ergs' 1 cm' 2 and for which accurate abun- 
dances could be derived. They have also excluded all galaxies with 
both [OIII]/l4959/H/3 < 0.7 and [OII],13727/H/3 > 1.0. Applying 
all these selection criteria, they obtain a sample of ~ 310 SDSS 
objects. Data for these sample objects include the emission line in- 
tensities of: [OIII]/14959,5007 A and [NII]^6584 relative to H/3 and 
the equivalent width of H/3. They also include the intensity of the 
[Oil] AA 3727,29 A emission line for the lowest redshift objects. 

The very young stellar population during the first Ma after 
each burst produces a high ionization, so the [OIII]/H/3 ratio is high, 
then it decreases raising again at 4 Ma due to the presence of Wolf 
Rayet s tars th at produces a harder continuum ^Garcia- Vargas et all 
Il995alll998h . After 5 Ma, the burst evolves to lower values of 
log[OIII]/H/3. Diagonal lines in this plot correspond to nearly con- 
stant ionization parameter and are swept along by the different suc- 
cessive bursts of a given age. High excitation objects, with high 
log[OIII]/H/3 and low log[OII]/H/3 ratios, have low metallicities 
and they are not reproduced by our model, which fits better the high 
metallicity sample, in the middle to right side of the panel. 

It can be better seen in Fig. [9] which shows the relation 
between the excitation parameter, log([OIII]/H/3), and the metal- 
licity indicator log([NII]/Ha) which gives infor mation about the 
young er populations plus the star formation rate dKennicutt et all 
1 19941) . At the beginning of the evolution, as we saw in previous 
sections, the star formation rate is high, and the Ho- emission of our 
model is strong. As the galaxy evolves, the Ha emission decreases 
while the [Nil] emission increases due to the growth of metallicity. 

Our models reproduce the general trend shown by data but 
cannot reproduce the lower metallicity (low [NII]/H„ ratio) objects 
to the left of the graph. This is due to the fact that the efficiency of 
the first burst, 33.1 per cent, produces a metallicity which is already 
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Figure 7. Time evolution of colours for model A2: a) U-B; b) B-V; c) V-R; d) V-I. The meaning of the (coloured) lines for different bursts is the same as in 
Fig. [5] Solid lines show the pure continuum colours while dotted lines also include the contribution by the strong emission lines. 




Figure 8. The relation between the oxygen line ratios for model A2. Solid coloured lines join the values reached at the same age by each different burst. Red, 
green, brown, blue and yellow colours correspond to ages 0-1, 2, 3, 4 and 5 Ma respectively.Different bursts, except the first one, are represented by different 
symbols: solid circles, squares, diamonds, triangles up, triangles left, triangles down and triangles right correspond to bursts at 1.3, 2.6, 3.9, 5.2, 6.5, 7.8 and 
9.1 Ga respectively. Open circles, squares and diamonds correspond to bursts at 10.4, 11,7 and 13.2 Ga respectively. The data are from references cited in the 
text. 



10 Martin-Manjon et al. 

1.5 1- 




I i i i i i l i l i ▼ l i l 

-3 -2.5 -2 -1.5 -1 -0.5 

log ([NH]/Hot) 

Figure 9. The relation between the excitation parameter, log([OIII]/H/3), and the metallicity indicator log([NII]/Ha) for model A2. The meaning of colours, 
symbols and lines is the same as in the previous figure. 



too high. A lower efficiency is therefore required to reproduce the 
observations of the less metallic objects. 

3.4 Combined continuum and emission line colours 

Once the continuum spectral energy distributions and the different 
emission line fluxes are computed, we can calculate the contribu- 
tion of these emission lines to the different broad band filters and 
then synthesize the colours of our model galaxies. These are the 
colours which are readily observable through integrated photome- 
try. 

We have done so taking into account only the strongest emis- 
sion lines that contribute to the colour in each broad-band spectral 
interval at redshift zero. These are: [Oil] AA 3727, 3729 A in U, H/5 
A 4861 A in B, [OIII]/li 4959, 5007 A in V, Ha A 6563 A in R and 
[SHI] AA 9069, 9532 A in I. 

The continuum colours and the colours including the emission 
line contribution at redshift zero are given in columns 3 to 10 of 
Table 1 for each of the models characterized by their burst number 
and age (columns 1 and 2). Columns 11 to 15 of the same table 
give the contribution, in percentage, by the emission lines to the 
different continuum fluxes. 

The necessary information to calculate these contributions for 
any other redshifts in the nearby universe (z <0. 1), is given in Table 
2 which lists for each of the models characterized by their burst 
number and age (columns 1 and 2), the continuum fluxes in each of 
the U,B,V,R and I (columns 3 to 7) and the fluxes in the [Oil], Hj8, 
[OIII], Ha and [SIII] lines (columns 8 to 12). 

Figure [7]shows, with dotted lines, the new computed colours 
including the contribution by the different emission lines. As it can 
be seen the colours do not change by a large amount in the case 
of U-B or V-R. However the effect in the B-V and V-I colours is 
rather dramatic. Figure [TOl shows the B-V vs V-I colour-colour dia- 
gram for model A2. Light blue circles show the continuum colours 
(solid lines in Figure |7]( while black solid ones show the colours 
computed including the contribution by the stronger emission lines 
given in Table 1. The inclusion of the contribution by the emis- 
sion lines to the continuum colours shifts the position of the model 



points almost perpendicularly to the originally computed ones. The 
location of blue points is mainly determined, even at zero redshift, 
by the contribution of strong [OIII]4959,5007 emissio n lines to the 
V band . In the same diagram we show the data by ICairos et al .1 
d2001alTbh on BCG as red squares with error bars. These correspond 
to readily observed colours therefore including both continuum and 
line emission. No reddening correction has been applied. It can be 
seen that some of the data are impossible to reproduce by the mod- 
els which do not take into account the contribution by emission 
lines, even if some amount of reddening is invoked (the arrow in 
the diagram shows the correction corresponding to 1 mg visual ex- 
tinction), while a certain amount of reddening would bring most of 
the data in agre e ment with the models. On he other hand the data 
by ICairos et all d2002l) of resolved locations in Mrk 370, shown 
in the figure as green triangles, correspond to data that have been 
corrected for the contribution by emission lines and therefore cor- 
respond to pure continuum colours. In these case, the agreement 
with the models is excellent. 



4 DISCUSSION 

The evolution of Hn galaxies must be considered by looking simul- 
taneously at both continuum and emission line properties, the first 
one providing information on a long time scale evolution, of the 
order of Ga, and the second one being related to the evolution in a 
muc h shorter time s cale, of the order of Ma. 

iDottoril dl98ll) suggested the use of the equivalent width of 
H/3, EW(H/?), as an age estimator for Hn regions, and applied this 
method to rank the ages of Hn regions in the Magellanic Clouds. If 
a SSP is considered, the H/? emission is very high at the beginning 
of the burst, while the continuum at H/3, dominated by the most 
massive and luminous stars, is low hence producing a large value 
of EW(H/3). As the burst evolves, the H/3 luminosity decreases as 
does the number of ionizing photons, and the contribution to the 
continuum increases thus lowering the value of EW(H/?). This case 
corresponds to the thick (black) line in Fig. II II 

The evolution of the equivalent width in our model is shown 
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Figure 10. B-V vs V-I colour-colour diagram for model A2. Light blue circles show the continuu m colours (solid line s in Figure |7J while black solid o nes 
show the colours computed including the contribution by the stronger emission lines. The data bv lCairos et"aT] feOOl j) on BCG and lCairos et alj J2002I) in 
Mrk 370 are included as red squares and green triangles respectively. 
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Figure 11. Time evolution of the H/3 equivalent width of model A2. The 
thick line corresponds to the first star burst of our model, a SSP with very 
low metallicity. Thin lines represent the successive bursts with same colours 
as in Fig. [5] 
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Figure 12. EW(H/3) as a function of the metallicity 12 + log(0/H) for the 
model A2. Colours have the same meaning than in Fig. [7] The first burst 
(solid black line) does not appeal' in the graph because its abundance is still 
the primordial one, 12 + log(0/H) « 1. 



in Fig.QT] It is seen that the first burst, that corresponds to a SSP, 
has an initial equivalent width slightly greater than the successive 
bursts, which decreases with time reaching very low values (< 40 
A ) 6 Ma after the beginning. When a new burst takes place, new 
massive ionizing stars form and the EW(H/?) goes up again, al- 
though not as much as in the previous one due to the decreasing 
gas mass involved in the burst, increased metallicity and the accu- 
mulated contribution from the continua from the previous bursts. 
At any rate, even taking into account this contribution, the EW(H/?) 
remains a good age indicator for the age of the current burst, inside 
10 Ma. On the other hand it is not possible to detect the presence of 
the underlying stellar populations for a given galaxy solely on the 
basis of this parameter. 

The distribution of the EW(H/J) in Hn galaxies shows that 



most of them have va lues lower than 150 A dTerlevich et al .LlT99ll : 
iHovos & DfazL l2006t) which, in our model s, corresponds to ages 
greater than 3-4 Ma. iTerlevich et alj | |2004|) . from the application 
of inversion methods, have shown that, globally, this distribution is 
inconsistent with what is expected from a single burst scenario and 
resembles more the distribution obtained for a succession of short 
bursts separated by quiescent periods with little or no star forma- 
tion. 

These results are not that strong when samples showing a re- 
stricted metallicity range are analyzed. However, as we have seen 
before, metallicity increases with time in a way which constrains 
the SFH which is best reproduced by our A2 model. Fig.ll2lshows 
the evolution of the EW(H/?) with oxygen abundance. During the 
duration of the first burst the EW(H/?) decreases while the oxygen 
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abundance increases from zero to 12 + log(0/H) ~ 8. The suc- 
cessive bursts start with this higher abundance and each of them 
evolves vertically in a 10 Ma time scale with the EW(H/3) decreas- 
ing rapidly at almost constant oxygen abundance. This occurs be- 
cause the oxygen is produced by the most massive stars in a very 
short time and therefore its abundance changes very little during the 
next 10 Ma after each burst (see Fig. [3}- It is clear that our model 
is ab le to reproduce t he mo st m etallic galaxies in th e samples taken 
from lHovos & Diaz! j2006l) and llzotov et"ai] d2006h where this high 
abundance is explained as the consequence of the gas being pro- 
cessed by one or more previous generations of stars. An alterna- 
tive scenario to produce such high oxygen abundances would be 
the occurrence of a very intense initial star formation burst, with 
a SFR of the order of 100 M Q Ma -1 , out of the range observed in 
local Hn galaxies (see Figures|3]and[T). On the other hand, galaxies 
12 + log(0/H) ^ 8 are not reproduced by our model. In order to 
fit these data, a model with a star formation efficiency lower than 
assumed, that is, < 33 per cent, should be used. 

The most informative data to uncover the presence of under- 
lying populations in Hn galaxies consist of the combination of a 
line emission parameter, which characterizes the properties of the 
current burst of star formation, and a continuum colour which rep- 
resents better the SFH in a longer time scale. Fig. Q/3] left panel, 
shows the equivalent width of H/3 vs the U-V colour along the evo- 
lution of the galaxy for the successive stellar bursts. When a given 
burst occurs the EW(H/3) is high and decreases as the stellar popu- 
lation ages. We hav e over-plotted the resu lts for the SSPs from the 
model Starburst99 dLeitherer et allll999l SB99) with two different 
metallicities as labelled. As expected, our results agree better with 
the higher metallicity SB99 model. The black line in the left panel, 
that represents a metal-poor SSP, is too blue compared with the 
data for a given EW(H/3). In order to decrease EW(H/?) and move 
the model colour to the red, a more metal-rich SSP was selected 
(dotted line at the left panel), but even such unrealistic high abun- 
dance does not reproduce the observations, as we have shown in 
Fig. [3] Therefore, the observed trend can not be explained as an 
age effect since EW(H/?) for low metallicity SSP does not reach 
values lower than 100A and does not reproduces the colors, nei- 
ther as a metallicity effect because such high abundances are not 
observed in Hn galaxies. It is therefore necessary to include both 
effects simultaneously, as we have done, to see how the effect of 
an underlying older population contributing to the color of the con- 
tinuum makes it redder. From the observational point of view, U-V 
colours for Hn galaxies are scarce, therefore w e have compared our 
results wit h those from | Hovos & Dfazl ( 2006) that i nclude galaxies 
taken from lTerlevich et alj Jl99ll) and lSalzer et al.l \ 19951) provid- 
ing the equivalent widths of both the [Oil] and [OIII] lines. We have 
then computed pseudo-colours from the intensities of the adjacent 
continua of [0II]/13727 and [OIII]^5007 lines as: 
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The right panel of Fig. Q/J] shows the equivalent width of H/3 
as a function of the pseudo-colour defined above for model A2 to- 
gether with the observational data. It should be noted that no cor- 
rection for extinction associated with the continuum light has been 



made. This correction amounts to about 0.15 dex for a standard 
reddening law and therefore a SSP would always be unable to re- 
produce the observed colours. Yet our Model A2 produce UV con- 
tinuum colours which are too blue as compared with the bulk of 
the observations which indicates that our assumed degree of atten- 
uation is too low. 



5 SUMMARY AN CONCLUSIONS 

We have explored the viability of a kind of a theoretical model in 
order to understand how the star formation takes place in Hn galax- 
ies. Our work involves the combination of three different tools: 
a chemical evolution code, an evolutionary synthesis code and a 
photo-ionization code, in a self-consistent way, i. e. taking the same 
assumptions about stellar evolution and nucleosynthesis, and takes 
into account the resulting metallicity in every time step. The chem- 
ical evolution code provides the star formation history of the galax- 
ies as well as the time evolution of the abundances in their interstel- 
lar medium, that can be confronted with observations. The evolu- 
tionary synthesis code provides the spectrophotometric evolution of 
the integrated stellar population. Finally, the photo-ionization code 
uses the properties of the ionizing stellar population to calculated 
the emission line spectra of the ionized gas. 

Each galaxy is modelled assuming an initial amount of unpro- 
cessed gas of 10 s M in a 1 kpc diameter. The evolution is com- 
puted along a total duration of 13.2 Ga during which 1 1 successive 
starbursts, separated by 1.3 Ga, take place . Different types of mod- 
els have been computed with different starbursting properties: equal 
and attenuated bursts, and different values of the initial star forma- 
tion efficiency. The only model which has proved to be able to fit 
the observational data is model A2, which consists of attenuated 
bursts with an initial star formation efficiency of 33. 1 per cent. This 
type of model reproduces the star formation rate estimates of local 
Hn galaxies and the oxygen abundances of the average metallicity 
objects although fails to account for the most metal deficient ones 
which would require models with lower star formation efficiencies 
than assumed. 

We have computed the time evolution of the continuum 
colours with and without the inclusion of the emission lines. It 
is shown that both nebular continuum and line emission must be 
taken into account in the photometric studies of the underlying stel- 
lar populations of Hn galaxies. The effect of the emission lines are 
more pronounced in the B-V and V-I colours. 

The combination of parameters which characterize the current 
star formation on a time scale of Ma, such as the equivalent width of 
H/3, and the SFH over a time scale of Ga, such as broad-band con- 
tinuum colours, is found to provide an effective means to uncover 
the presence of underlying stellar populations. The comparison of 
models and observations show that in most Hn galaxies SSP are 
unable to reproduce the relatively red colours shown by data with 
the observed H/? equivalent width values requiring the contribu- 
tion of previous stellar generations. Our model, however, produces 
U-V colours that are too blue when compared with observations, 
therefore implying that the successive bursts that take place in Hn 
galaxies should be even more attenuated than has been assumed in 
our work. 
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Figu re 13. Lef t panel : The equivalent width of Hf! as a function of the continuum colour for model A2. The results for the SSPs from the model Starburst99 
iLeitherer et al ] fl999h with two different metallicities as labelled are overplotted. Colours have the same me aning as in Fig. \7\ Rig ht p anel: The equivalen t 
width of H/3 as a function of the pseudo-colour defined by equations 13 to 15 for model A2. The data are from lTerlevich~ et al. 1 1991) and lSalzer et"aU jl995h . 
In both panels the large dots mark ages separated 2 Ma. 
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Tab] 

the 1 



le 1. Evolution of continuum and total broad band colors for different bursts along the time after each stellar burst occurs. Only 
last burst evolution is shown. The whole table is available in electronic format. 



NB 





(U-B) c 


(B-V) r 


(V-R)c 


(V-I) c 


(U-B) 


(B-V) 


(V-R) 


(V-I) 


%U 


% B 


% V 


%R 


% I 


5.75 


-1.205 


0.019 


0.271 


0.480 


-1.292 


0.748 


0.109 


0.064 


21.254 


14.681 


57.230 


50.400 


37.264 


5.80 


-1.202 


0.015 


0.267 


0.474 


-1.300 


0.694 


0.147 


0.102 


21.845 


14.497 


54.820 


49.812 


36.728 


5.85 


-1.203 


0.014 


0.267 


0.472 


-1.301 


0.706 


0.142 


-0.065 


21.900 


14.506 


55.610 


50.230 


42.959 


5.90 


-1.203 


0.012 


0.264 


0.468 


-1.311 


0.678 


0.164 


0.117 


22.634 


14.424 


54.460 


50.115 


37.121 


5.95 


-1.200 


0.008 


0.260 


0.461 


-1.345 


0.611 


0.222 


0.177 


25.040 


14.260 


51.674 


50.000 


37.274 


6.00 


-1.198 


0.004 


0.255 


0.453 


-1.378 


0.555 


0.269 


0.225 


27.207 


13.997 


49.152 


49.839 


37.295 


6.05 


-1.197 


0.001 


0.252 


0.448 


-1.383 


0.535 


0.281 


0.236 


27.391 


13.889 


48.246 


49.628 


37.110 


6.10 


-1.196 


-0.002 


0.248 


0.441 


-1.397 


0.500 


0.307 


0.260 


28.277 


13.716 


46.603 


49.456 


36.942 


6.15 


-1.191 


-0.009 


0.238 


0.426 


-1.435 


0.407 


0.374 


0.327 


30.679 


13.226 


41.893 


48.686 


36.373 


6.20 


-1.189 


-0.015 


0.232 


0.414 


-1.457 


0.355 


0.412 


0.363 


31.957 


12.892 


39.095 


48.422 


36.106 


6.25 


-1.185 


-0.021 


0.224 


0.401 


-1.477 


0.300 


0.452 


0.399 


33.171 


12.523 


36.012 


48.132 


35.917 


6.30 


-1.172 


-0.035 


0.202 


0.367 


-1.517 


0.136- 


0.554 


0.504 


35.548 


11.422 


25.592 


46.160 


34.394 


6.35 


-1.160 


-0.047 


0.182 


0.335 


-1.515 


0.027- 


0.601 


0.551 


35.367 


10.370 


17.718 


44.104 


32.609 


6.40 


-1.138 


-0.067 


0.148 


0.278 


-1.475 


-0.068- 


0.608 


0.551 


33.098 


8.712 


9.986 


41.101 


30.006 


6.45 


-1.113 


-0.083 


0.116 


0.226 


-1.403 


-0.107- 


0.543 


0.483 


29.051 


7.251 


6.539 


36.974 


26.347 


6.50 


-1.073 


-0.107 


0.061 


0.128 


-1.274 


-0.152- 


0.407 


0.336 


20.862 


4.735 


2.161 


28.889 


19.217 


6.60 


-1.020 


-0.076 


0.084 


0.182 


-1.079 


0.316- 


0.000 


0.041 


12.736 


4.400 


32.054 


24.293 


16.869 


6.70 


-0.878 


-0.006 


0.119 


0.277 


-0.976 


0.024- 


0.196 


0.272 


10.700 


2.135 


6.592 


13.066 


6.218 


6.80 


-0.484 


0.174 


0.206 


0.446 


-0.497 


0.140 


0.238 


0.447 


1.707 


0.518 


0.000 


2.901 


0.239 


6.90 


-0.553 


0.244 


0.291 


0.617 


-0.553 


0.212 


0.304 


0.617 


0.295 


0.265 


0.000 


1.341 


0.018 


7.00 


-0.630 


0.219 


0.331 


0.723 


-0.628 


0.190 


0.339 


0.723 


0.029 


0.146 


0.000 


0.706 


0.001 


7.10 


-0.567 


0.210 


0.152 


0.698 


-0.567 


0.182 


0.314 


0.698 


0.000 


0.059 


0.00 


0.00 


0.00 



Table 2. Evolution of broad band and the strongest emission lines luminosities along the time after each stellar burst occurs. Only 
the last burst evolution is shown. The whole table is available in electronic format. 



NB log | L(/ L B L v h R U L([OII]) L(H/3) L([OIII]) L(Ha) L([SIII]) 

10 40 erg s -1 



11 5.75 


1.559E+01 


1.423E+01 


7.080E+00 


1 


123E+01 


7.442E+00 


4.252E+00 


2.474E+00 


9.574E+00 


1.1533E+01 


4.466E+00 


11 5.80 


1.569E+01 


1.435E+01 


7.126E+00 


1 


126E+01 


7.442E+00 


4.431E+00 


2.458E+00 


8.812E+40 


1.1296E+01 


4.364E+00 


11 5.85 


1.582E+01 


1.447E+01 


7.172E+00 


1 


134E+01 


7.483E+00 


4.483E+00 


2.480E+00 


9.O8OE+O0 


1.1560E+01 


4.459E+00 


11 5.90 


1.595E+01 


1.460E+01 


7.225E+00 


1 


139E+01 


7.510E+00 


4.716E+00 


2.487E+00 


8.731E+00 


1.1560E+01 


4.480E+00 


11 5.95 


1.612E+01 


1.478E+01 


7.286E+00 


1 


144E+01 


7.524E+00 


5.439E+00 


2.483E+00 


7.872E+00 


1.1560E+01 


4.518E+00 


11 6.00 


1.631E+01 


1.500E+01 


7.367E+00 


1 


151E+01 


7.552E+00 


6.160E+00 


2.466E+00 


7.196E+00 


1.1560E+01 


4.538E+00 


11 6.05 


1.654E+01 


1.519E+01 


7.448E+00 


1 


161E+01 


7.601E+00 


6.303E+00 


2.476E+00 


7.016E+00 


1.1560E+01 


4.532E+00 


11 6.10 


1.680E+01 


1.545E+01 


7.545E+00 


1 


172E+01 


7.650E+00 


6.692E+00 


2.481E+00 


6.654E+00 


1.1586E+01 


4.528E+00 


11 6.15 


1.714E+01 


1.584E+01 


7.685E+00 


1 


184E+01 


7.685E+00 


7.666E+00 


2.439E+00 


5.598E+00 


1.1349E+01 


4.439E+00 


11 6.20 


1.751E+01 


1.622E+01 


7.835E+00 


1 


199E+01 


7.756E+00 


8.309E+00 


2.425E+00 


5.082E+00 


1.1375E+01 


4.429E+00 


11 6.25 


1.792E+01 


1.666E+01 


8.OO3E+O0 


1 


216E+01 


7.821E+00 


8.986E+00 


2.410E+00 


4.551E+00 


1.1401E+01 


4.429E+00 


11 6.30 


1.856E+01 


1.746E+01 


8.281E+00 


1 


234E+01 


7.850E+00 


1.034E+01 


2.275E+00 


2.878E+00 


1.0689E+01 


4.158E+00 


11 6.35 


1.927E+01 


1.833E+01 


8.599E+00 


1 


257E+01 


7.908E+00 


1.065E+01 


2.143E+00 


1.871E+00 


1.0021E+01 


3.866E+00 


11 6.40 


2.050E+01 


1.990E+01 


9.172E+00 


1 


299E+01 


8.OO3E+0O 


1.024E+01 


1.919E+00 


1.028E+00 


9.1614E+00 


3.467E+00 


11 6.45 


2.172E+01 


2.158E+01 


9.810E+00 


1 


349E+01 


8.152E+00 


8.986E+00 


1.704E+00 


6.936E-01 


7.9977E+00 


2.946E+00 


11 6.50 


2.478E+01 


2.554E+01 


1.136E+01 


1 


485E+01 


8.631E+00 


6.600E+00 


1.282E+00 


2.535E-01 


6.0948E+00 


2.074E+00 


11 6.60 


1.929E+01 


2.088E+01 


9.543E+00 


1 


273E+01 


7.615E+00 


2.165E+00 


7.670E-01 


3.724E+00 


3.4910E+00 


1.375E+00 


11 6.70 


1.260E+01 


1.556E+01 


7.566E+00 


1 


043E+01 


6.596E+00 


1.526E+00 


3.431E-01 


5.396E-01 


1.5847E+00 


4.419E-01 


11 6.80 


1.012E+01 


1.795E+01 


1.023E+01 


1 


528E+01 


1.041E+01 


1.776E-01 


9.445E-02 


0.000E+00 


4.6128E-01 


2.525E-02 


11 6.90 


7.074E+00 


1.177E+01 


7.146E+00 


1 


154E+01 


8.520E+00 


2.121E-02 


3.171E-02 


0.000E+00 


1.5847E-01 


1.615E-03 


1 1 7.00 


5.O45E+O0 


7.821E+00 


4.648E+00 


7 


799E+00 


6.110E+00 


1.491E-03 


1.159E-02 


O.OOOE+00 


5.6100E-02 


8.127E-05 


11 7.10 


4.390E+00 


7.205E+00 


4.247E+00 


7 


009E+00 


5.456E+00 


6.406E-08 


4.361E-03 


0.000E+00 


2.089 1E-05 


0.000E+00 



